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ABSTRACT 



Three regression techniques for the calculation of ballistic winds 
using satellite infrared spectrophotometer (SIRS) data are evaluated. 
Experiments testing the derived methods are carried out using two 
SIRS-B data set. In all three cases, the data sets are divided 
into a large dependent sample used to derive the regression equations 
and smaller independent sample for testing the equations' validity. 

One method, that of geostrophic calculations using regression 
derived "D" values of ballistic geopotential height is found to be 
superior to the other methods. 
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I. INTRODUCTION 



The resultant wind which deflects a reentry vehicle as it 
travels along its flight path is known as the ballistic wind. This 
is a weighted integral of the winds encountered at the various levels 
of the atmosphere through which the projectile travels. The general 
form of the ballistic wind calculation is: 




, ^0 n p) W, (l n p) c!0n p) 



0 ) 



where is the ballistic wind, p is pressure, p qq is surface pressure, 
W is the actual wind at level p, and W is the ballistic wind 
weighting factor for a layer. W is proportional to the effect of 
the winds in a particular layer and the time the vehicle remains in 
that layer. A similarly computed quantity, ballistic density, is 
representative of the total retarding force of the atmosphere on a 
reentry vehicle. 

Present methods of evaluation of ballistic winds and density 
involve the use of rawinsonde data and rocketsonde data as available 
for the levels in the integration of Eq. (l). Smith, et. al. [6] 
have recently reported success in obtaining atmospheric temperature- 
pressure profiles using satellite infrared spectrophotometer (SIRS) 
data. This success implies that ballistic winds and densities could 
be calculated in any region without the necessity of a rawinsonde 
release. 
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Ruggles [5] has noted that the envelope of weighting functions 
used with the individual SIRS channel radiances in obtaining 
atmospheric soundings and the ballistic wind weighting factors are 
analogous. He further suggested that the ballistic winds and 
densities might be calculated directly from SIRS data without first 
calculating the temperature-pressure profile. Ruggles' development 
is outlined below. 

It is the object of this research to investigate the feasibility 
of Ruggles' suggestion with respect to ballistic winds. 
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II. DEVELOPMENT OF EQUATIONS 



A. FEASIBILITY OF THE USE OF SIRS DATA 

If the geostrophic approximation is assumed, Eq, (l) can then 
"be expressed as: 




--f-IK *%Z(Vj ) '»p)V/ v (lnp) d(Up) 

o J 



( 2 ) 



or, rearranging terms: 



k ik * V|Q r* **■ { * p) h (^i°) ^ ipi (p) 



where g is the acceleration due to gravity, f is the Coriolis 

parameter, andtf^ is the horizontal gradient of geopotential height 

on a constant pressure surface. Now, expressing the weighting 

function W as a linear combination of known arbitrary functions 
v 

of In p, F^ (in p) : 



n 



W y (In p) = £ a i F ± (in p) 



(4) 



1 = 1 



where are constants , Similarly, Z(x, y, In p) can be expressed 
as: 



n 

Z(x, y, In p) = ^ r ± H ± (x, y, In p) (5) 

i=1 

where r^ are the constants and H^ are a set of arbitrary functions 
of In p. Substitution of Eq. (4) and Eq. (5) into Eq. (3) produces: 
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( 6 ) 




-9 

f Inp 

r»o 




Now consider the radiative transfer equation for the radiance 
observed by the i^* 1 channel of the SIRS, sensing at a wave number, l) . 

w ; = Wo [^V J TCKp\X(VJ^]4 0 B.&) ) TM J(l^ (7) 

P.o ei(in^) 

Here, N. is the radiance sensed, and N,. is the radiance emitted 
’ l ^i 

by a contaminant such as the top surface of clouds. This N n is a 

u i 

function of A, the amount of cloud cover, , the wave number of the 
i^* 1 channel, and T(ln p) the temperature structure above the clouds. 
The Planck radiance at wave number, V , and temperature, T, is 
represented by B.. LV , T(1 n pj) , an bTCP, In p) is the fractional 
transmittance of the atmosphere at wave number from pressure level 
p to the SIRS. 

Smith et. al. [6] have shown that cloud contaminant term can be 
evaluated by an iterative procedure. With that term evaluated, 

Eq, (7) becomes: 

K *i -I MV,TO*f0 g j fe k.r) d(| np ) (8) 

Po ° d (l*p) ' 

•X- 

where N . is the radiance corrected for cloud contamination effects, 

l 

also called the clear column radiance. Since the only constraints 
on F^ are that they are known functions of In p, we can specify, 

(after Ruggles [ 5 ] ): 
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dT 1 (V, ^ p) 

p = = 

i d In p 

Similarly, since B.[V> T(ln pj is an increasing function of T and 
also of In p, it satisfied the constraints put on IL in Eq. ( 5 ) and 
we can set: 



H. = B.[V, T(ln p)] 

Substituting these expressions into Eq. (8) produces: 

h'fhMCM (?) 

Using this in equation (6) results in a direct expression for ballistic 
winds in terms of the SIRS radiances. 
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IK x s? ( a. r N. 

p <s ^ 



I 



( 10 ) 



Now n is the number of SIRS channels used. In practice, the constant 
terms are grouped together and are statistically determined through 
the use of the BIHED 02R stepwise multiple regression routine, (Dixon, 
1 ) available at the ¥. R, Church Computer Center at the Naval 
Postgraduate School. This produces a final equation of the form: 



K="f lkxv p 





( 11 ) 



where b^ are the regression coefficients and is the regression 
equation constant. The constant terms, g and In are absorbed by 
the coefficients, b., 

7 i 
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B. "D" VALUE APPROACH 



The bracketed term in Eq. (ll) represents the regression equation. 



If the ballistic wind, V , is geostrophic, this term can be looked 

JJ 

upon as a ballistic geopotential height. The regression equation can 
be derived by correlation of SIRS radiances with a dependent sample 
of ballistic geopotential heights. These are obtained by computing 
a weighted average of observed geopotential heights in the vertical. 
The ballistic wind weighting factors are used in the averaging. For 
the dependent sample, deviations from standard, or "D" values of 
height were used for all available mandatory levels from rawinsonde 
soundings. For levels above the maximum height attained by the 
rawinsonde (usually 10 mb) temperature values were extrapolated using 
the standard lapse rate for 45N in January, The hydrostatic 
approximation was then -used to calculate ”I) ,! values. Linear inter- 
polation was used in the rare instances when an interior mandatory 
level value was missing in a sounding. The following equation 
illustrates the method used for obtaining the dependent sample of 
"D" values. 



Here. Ik, is the ballistic "L" value, ¥ are the ballistic wind weighting 



actual "D" values of height at each mandatory level. 

This method of calculating \V then, consists of two steps. 

JJ 

First a regression equation is derived by correlation of SIRS radiance 



\ rt \\ o 



d e = £ Wj, 



p p 



( 12 ) 
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data with the dependent sample of "D" values calculated using Eq. 
(12). Then, Eq. (ll) is used with the derived regression equation 
to calculate the ballistic wind. 



C. GRADIENT OP RADIANCES APPROACH 

Rearrangement of the order of Eq. (11) yields 



' V 4~ [ £ c ; IKxv « N ; +c z] 



(13) 






The regression coefficients are represented by c^, and C ^ is the 
regression equation constant. Here, c^ and C^ are different than 
b^ and C^ in Eq. (11). Though equivalent, this is procedurally 
different in evaluation from Eq. (11). In this case the entire 

expression is a regression equation. The equation is derived by 

I \rJ. 

correlation of values of -r- PIP- with v, the north-south comoonent 

$ d* B 

of V-g. Similarly, Ug, the east-west component is correlated with 

\ AN* 

# The dependent sample of vu. and v^ is obtained as before 

J dy b b 

from rawinsonde soundings using Eq, (1), Again, missing interior 
mandatory level values were supplied by linear interpolation. 
Missing upper level values of u were extrapolated by paralleling 
a mean wind profile. Missing v values were extrapolated by simply 
decreasing the already small v component by five per cent with each 
level increment. This was justifiable since v components above 
10 mb were small in all cases. 



D. DIRECT APPROACH 

A purely statistical method suggested by Dr, F, L, Martin was to 
perform direct correlations between the SIRS radiances and the ballistic 
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wind components. This produces a very simple regression equation 
for the ballistic winds. 



\ = ( W. N. + K ( 14 ) 

1^.1 



Here, hu are the regression coefficients and K is the regression 
equation constant. 
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III. SIRS-A EXPERIMENT 



A. DATA AND DATA HANDLING 

The data used in the SIRS-A experiment were obtained from Project 
FAMOS and consisted of radiances, corrected for cloud contamination, 
for the eight SIRS channels for two four-day periods: December 27-30, 

1969, and January 2-5, 1970. The area of data coverage was the 
European continent including the area of the U.S.S.R. west of the 
Urals, and the time period covered was from 0800 GMT to 1200 GMT for 
each day. The dependent ballistic "D" values and winds were computed 
from the 1200 GMT rawinsonde data from the same periods. The 
rawinsonde data did not include any data from France or the British 



satellite subtracks lying in this data sparse zone. In addition, data 
was typically quite sparse in the Mediterranean region. 

The 27 degree longitudinal spacing of Nimbus III SIRS data 
precludes the calculation of both u and v components of ballistic 
winds by the geostrophic approximation. Instead, Vg, T , the wind 
component normal to the satellite subtrack was calculated using 
Eqs, (11) and (13) modified as indicated below 



Isles. In some cases thrs was a Iruixtation due to the locate. 




(11a) 
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Here n is the generally northerly direction parallel to the satellite 
subtrack. The equation for the direct approach remained unchanged. 

The "D" values and values of u^ and v.g calculated from all 
available rawinsonde data were plotted and each of the fields analyzed. 
Values were then interpolated to SIRS subtrack positions. Since the 
points with SIRS data were irregularly spaced, the data were processed 
using a cubic spline fitting interpolation routine which yielded 
radiance data at evenly spaced points at a 2 degree latitude interval 
along the satellite subtrack. The effects of this interpolation on 
the methods’ accuracy are discussed later. A typical Nimbus III 
subtrack and data is shown in Pig. 1 . 

B. "D” VALUE APPROACH 

The initial regression equations used January 2-4 data as a 
dependent data set. Initial experimentation was confined to the 
January data because the missing French and British data caused a 
bad gap in data coverage relative to the positions of SIRS data on 
December 29* 

At first, two dependent samples of "D" values were correlated 
with radiance data, one using interpolated evenly spaced radiance 
data and another using data from the uninterpolated sub-points. 

Results using the two different data sets were similar and are 
summarized in Table 1 . Since the difference in results was slight 
and the evenly spaced data were much easier to handle, experimentation 
proceeded with the interpolated data. 
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© UNINTERPOLATED SIRS-A 
DATA POINTS 

© INTERPOLATED DATA POINTS 

V Dld COMPONENTS 

B N 



Fig. 1 . Typical Nimbus III Satellite Subtrack Showing 
Uninterpolated ana Interpolated SIRS-A Data. 
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TABLE I 



Effect 


of Interpolation < 


on Results of "D" 


Value 


Correlations with 


SIRS Radiances 




Correlation 


Standard Error of 




Coefficient 


"D" Value Estimate 


Interpolated 


Data 


.9 75 


51.6 meters 


Unint e rpo lat e d 


Data 


.979 


49.7 



Since the standard errors of "D" value determination were quite 
large compared to those found by Elsberry and Martin [ 2] for a 
dependent sample of simulated radiance data (22.3 meters), the data 
were carefully examined. It was noted that the great majority of 
points which nad large residuals between rawinsonde derived "If ; 
values and those obtained from the regression equation were located 
in data-poor regions. It was therefore assumed that the large 
residuals were mainly due to inaccuracies in the dependent "D" 
values. 

In an attempt to reduce the standard errors of the estimated 
"D" values, those data points which had produced residuals greater 
than 50 meters were removed from the sample. For this filtered 
sample, a new correlation was calculated with somewhat improved 
results, but this reduced the size of the data set from n = 90 to 
n = 49 • The correlation was increased to .987* and the standard 
error was reduced to 28.3 meters. 

For the December data, the same two correlation runs were 
performed with the following results: for all data points (n = 80) 
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the multiple correlation coefficient was .806 and the standard 
error was 51*8 meters. For the filtered data (n = 53) the correlation 
was .946 and the standard error was 27 .6 meters. The four regression 
equations for ballistic "D" value, Dg, are listed below: 

December, unfiltered data: 



D^. = 1.281N, - 5 . 243 N_ + 1 5 . 642N + 7.884N. - 58.512N C 
r> \ 5 4 3 

+ 63.779N 6 + 9.836N - 1.088N Q - 1822.8 



(15) 



December, filtered data: 

D- = 1.2121L - 7-702N o + 14.870N 2 + 8.320IT - 58.870N C 

D 1 2 5 4 3 

+ 66.715N 6 + 12 . - 1.834N 8 - 1986.0 



( 16 ) 



January, unfiltered data: 



D b = -5.497N 1 - 58.457N 2 - 42 . 4841*3 + 52.944^ - 91.296^ 
+ 40.878N 6 - 10.980N 7 + 14.980N q + 1792.3 



( 17 ) 



January, filtered data: 



d b = -5.622^ - 51*070N 2 + 28.321N 5 + 20.475^ - 50.746^ 

+ 19.161N 6 + 1.928N 7 + 12.966N q + 1275.7 



(18) 



A comparison of the magnitudes of the regression coefficients 
indicated the significance of the different SIRS channels varies 
markedly from the December to the January data sets. In December, 
channels 5 and. 6 provide the major modification to the constant 
term while in the January case, channels 2, 3> 4> 5> and 6 all had 
contributions of roughly the same order of magnitude. The coefficient 
for channel 5> however, is almost twice as large as the coefficients 
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for the other channels. Results reported by Smith, et. al. [6] would 
indicate that the principal contribution in both cases came from the 
layer centered around 25 O-I 5 O mb. In the January case, however, it 
would appear the principal contribution came from a slightly lower 
level and that the vertical wind shear was less. 

That the regression equations derived for the filtered and 
unfiltered cases were nearly equivalent was demonstrated by comparing 
"D" values calculated from both periods. The comparisons are 
tabulated in Table II. 



TABLE II 

Comparison of "D" Values Obtained from Filtered 
and Unfiltered Regressions 





Correlation 


Explained Variance 


Sta. error 


Dec. 


.993 


98.6$ 


8.7 meters 


Jan. 


.995 


99.0 


20.5 



To check this scheme for calculating usable ballistic winds, 
the "D" values calculated using all four of the derived regression 
equations, Eqs. ( 15 ) - (l 8) were used with Eq. (lib) to generate 
ballistic wind values. These were then compared with the rawinsonde 
derived values as summarized in Table III. 

The mean of the absolute values, JV^|, of the wind components is 
included in the table since both easterlies and westerlies may be 
included in the sample. In December, the winds were quite variable 
and W-g,, and (\, r _j,| differ greatly. In January the winds were highly 
zonal and andfVg^j are equal. 
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Here, VI represents ballistic winds calculated using the regression 
equation derived from all data points, and V2, those obtained from 
the equation derived from filtered data. The (i) indicates the 
independent samples, either December 30 or January 5* Although the 
filtered and unfiltered dependent equations were used, both were 
applied to the same, unfiltered data set for the test of the independent 
sample. An unfiltered data set was used as a check since it would 
be impossible to perform the filtering process on an independent 
data set if the procedure were to be used operationally. 

It was found that the significance of the channels entering into 
the dependent regression equations, Eqs. ( 1 5 — 18) , became very slight 
after more than five channels were used. The correlation coefficient 
was increased by less than 1.0 per cent with the addition of each 
succeeding channel. In spite of this, regression equations were 
derived using all SIRS channels since the addition of succeeding 
channels reduced the standard error. 

It can be seen that the filtering of the data had little effect 
on the final outcome. In fact, the results were opposite in the 
two cases. For the independent December data, filtering slightly 
improved the results, while it slightly degraded them in January, 

In either case, the difference was insignificant, at most amounting 
to a 6.1 per cent difference in explained variance or a 2.3 knot 
difference in standard error. The standard errors of the wind 
component are not overly large considering the wind is a weighted 
average from the surface to 30 km. In any case, they are well 
within the accuracy criteria set down by Finke [ 3 ]* 
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For both periods, the wind estimates obtained by using the 
geostrophic approximation with SIRS derived "D" values are better 
than those derived from "D" values obtained from the rawinsondes. 
Geostrophic values of V ^ were calculated using "D" values interpolated 
to the satellite subpoints and compared with the calculated values 
of V-™ at those points. A correlation of .725 was obtained for the 
December data and .766 was the correlation for January. An explanation 
of the superior performance of the SIRS techniques is that the "D" 
value fields obtained from SIRS data are free of the inaccuracies 
inherent in a field obtained by many different sensors. Unlike 
separate rawinsondes, the SIRS instrument tends to make all errors 
in the same direction, thus not significanly affecting gradients. 

C. GRADIENTS OF IT, 

i. 

This approach is somewhat more direct than the "D" value approach. 
Here, the derivatives ~ , were compared directly with the 

5 dr\ 

calculated values of V- and the resulting regression equations were 

-D 

used to calculate V- for the independent sample. The regression 

JJ 

equation derived from the December data is given below as an example. 



V BN = " 2 '3T1x 1 - 11.456x 2 + 30.703x 3 - 1.548x 4 

+64.865x 5 - 25.090x 6 - 2.122x ? + 2.251x Q + 24.42 



(19) 



Here, x. indicates “ 7 - 1 . 

This approach was tried with data interpolated to evenly spaced 
grid points and also with the derivatives calculated at the mid-points 
between the actual, uninterpolated satellite subpoints. In view of 
the results obtained during the first phase of the experiment, no 
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data filtering was attempted. The results of this approach are summarized 
in Table IV, Tests of the December 30 and January 5 independent 
samples for the in situ SIRS data were not made in view of the more 
favorable results for the dependent data when subtrack interpolation 
was employed. 

The results here are somewhat mixed. Only for the December 30 
case does the scheme show as much accuracy as the results for the 
M D M value approach shown in Table III, The differences between the 
results of this method and those of the "D" value approach may be 

explained by the greater number of derivatives which were taken in 

1 

the -r 2 — c approach. In the latter case, a derivative was 

1 or\ 

calculated for each SIRS channel, while in the "D" value approach 
the only derivative calculated was that of the ballistic M D" value. 

The smaller explained variance with the ■“ approach probably 

resulted from a summation of the roundoff errors in each derivative 
calculation. In this approach, interpolation to obtain uniform 
spacing of data produced a greater improvement in results than it did 
in the "D" value approach. This is perhaps explainable by the 
very small values of compared to the magnitudes of , With 

such small values, the loss of definition due to the wide spacing of 
some of the uninterpolated data could seriously degrade results. The 
improvement resulting from interpolation is indicative of the cubic 
spline fitting routine’s ability to accurately determine the data 
patterns. 
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D. DIRECT APPROACH 



A third method tried was the derivation of a regression equation 
by direct correlation between Vg and 1-L . This approach gave poorer 
results than the previous two methods used, as Table V' indicates. 
(See also Tables III and 1?.) 

TABLE V 

Comparison of Winds Obtained Prom Direct 
Approach with Rawinsonde-Derived Winds (SIRS-A) 







Corr 


Expalined 


Standard 


Max 


Sample 


n 


Coeff 


Variance 


Error 


Error 


27-29 Dec 


80 


.788 


62.276 


12.8kts 


34.3kts 


30 Dec 


31 


.771 


59*5 


18.0 


31*1 


2-4 Jan 


88 


.607 


36.8 


15.6 


38.8 


5 Jan 


23 


.372 


13*8 


13*3 


21.1 


Except in 


the 


case of 


the December- 


27-29 "D" value 


: approach 


these results 


are 


somewhat 


inferior to 


those obtained by either 



the first two methods described. Two factors could account for the 
superior performance in December. First, the value of \Vg is smaller 
in December which would tend to make smaller errors. Secondly, 
the patterns obtained on the u and v analyses were more constant 
with time in the December case and thus the regression coefficients 
obtained from the dependent sample would be more applicable to the 
independent data. 
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E. SUMMARY 



Both the "D" value approach and the gradient of iL approach 
appear capable of providing usable values of \V_,, The direct 

-D 

approach, however, in most cases gave correlations too low to be 
considered of value. It is interesting to note that the correlations 
for the independent data sets were higher than for the dependent 
data with the December data. This result was also noted by Elsberry 
and Martin [ 2] in their use of the same data for ballistic density 
calculations. 

To further check the temporal stability of the techniques, the 
December equations were applied to the January data set and vice 
versa. The resulting values of ballistic wind were then compared to 
the rawinsonde-de rived values. Table VI gives a summary of the 
results for the ,! D" value approach. All runs were performed using 
using equally spaced data. 



TABLE VI 

Results of Temporal Stability Checks 
of the "D" Value Approach 



Sample 

Dec Data 
Jan Equations 


Corr 

Coeff 


Explained 

Variance 


Std 

Error 


VI -V 
BN 


.123 


1.5$ 


21 .Okts 


v r Y BN^ I),I)ec 50 


.677 


45.9 


22.0 


Jan Data 
Dec Equations 








VI -V 

BN 


.413 


17.1 


13.1 


V1 " V BN (l) ’ Jan 5 


.370 


13.7 


11.2 
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It can be seen that the temporal stability is rather poor. 

The percentages of explained variances were generally so low as 
to be insignificant. It is interesting to note that the best 
performance was for the December JO data using equations formed from 
the later data set, A possible implication to be drawn from Table 
VI is that the regression equations would have to be continually 
updated using data from a period immediately prior to the day for 
which the independent wind values are to be calculated. 

In the tables of results for the various methods, it may be 
noticed that the sample sizes vary considerably although the samples 
are taken from the same data set. This is the result of the loss 
of end points from the rows of data in the finite differencing 
processes used to calculate derivatives. 
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IV. SIRS-B EXPERIMENT 



A. DATA AND DATA HANDLING 

The SIRS-B data used in the second phase of experimentation were 
taken from the period January 10-17, 1971# Data were selected from 
the latitude belt from 40N to 75N. As was the case in previous 
experimentation, corrected clear column radiances were used. 
Concurrent rawinsonde data were also obtained. Though it was 
expected to sacrifice some data density to obtain the side-looking 
capability of the SIRS-B, data coverage vas still generally poor. 
Data points were widely scattered and many areas of good rawinsonde 
coverage were completely devoid of the SIRS-B data. In areas with 
SIRS— B coverage, the normal spacing between data points vas 450 to 
65 O kilometers. Another shortcoming of the SIRS-B data was that 
the window channel was inoperative so that only seven channels 
were available. 

To permit the calculation of gradients, the data were chosen 
from areas where several SIRS data points could be found grouped in 
a "data cluster." These clusters consisted of one master data 
point surrounded by three to six other "slave" points. A schematic 
of a data cluster showing a "master" and five "slaves" is shown 
in Fig. 2. The maximum distance between a master and a slave was 
arbitrarily limited to 1000 kilometers. The surface between each 
master-slave pair was assumed to be linear and a plane was fitted 
between each pair. The average slope for all of the master-slave 
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@ MASTER 
O SLAVES 



Fig. 2 Sample SIRS-B Data Cluster 
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pairs was then taken as the partial derivative with respect to 
the x or y direction for either the "D" value or the radiance 
value at the master point. 

For the period January 10-17, 45 separate data clusters were 
found in areas with good rawinsonde coverage. These clusters 
comprised a total of 194 data points. The first 30 clusters 
covering January 10-14 were used as the dependent sample in all 
experiments and the 15 clusters from January 15-17 were used as 
the independent sample. The side-looking capability of the SIRS-B 
instrument enabled calculation of both the u and v components of 
the ballistic wind. Eqs. (ll) and (13) were therefore broken 
down into their appropriate components. Otherwise, computational 
procedures used were essentially the same as those used in the 
SIRS-A experiment. 

B. THE "D" VALUE APPROACH 

The regression equation for calculation of ballistic "D" 
values was obtained by correlation of rawinsonde-derived "D" 
values with the radiance data for all 1 30 data points in the 
dependent sample. The resulting equation is given as Eq. (20). 

D b = -10.82N 2 + 15.23N 4 - 47.41N 6 + 95.96n ? - 30.45N Q 

- 1727.6 

Note that in this equation, channels 3 and- 5 were not statistically 
significant. 




The correlation coefficient between "D" values and radiances 
for the dependent sample was .$ 06 , indicating 81.2 per cent of the 
variance was explained. The standard error was 119*8 meters. Though 
this standard error was quite large, no data filtering was attempted 
due to the already small sample size and the negligible improvements 
with filtering in the SIRS-A experiment. 

"D" values obtained from the application of Eq. (20) to the data 
points in the indpendent sample were next used in Eq. (ll) to obtain 
u and v components of ballistic wind. These SIRS-3 derived ballistic 
winds were then correlated with rawinsonde ballistic winds to check 
the accuracy of this method. The results of these correlations for 
the independent sample are tabulated in Table VII. 

It can be seen that the "D" value approach does quite well at 
determining the v component of \V _. In fact, the results are 

-D 

comparable to those for the one-day January dependent sample as noted 
in Table III. The poorer results with the u component may have 
resulted from the relatively small population of sample clusters. 

The method used to calculate gradients, when used over the wide 
spacings involved in a typical cluster might also contribute to 
the low correlations. 



C. GRADIENTS OP IT 



This technique, using Eq. ( 1 3) » was carried out in essentially 



the same manner as the experiment with SIRS-A data. In this case, 

\ \ 

however, the u and v components were correlated with values of 

and -7* respectively. 

7 

The resulting equations for these components, u^ and v T are given 
below. 
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Comparison of Ballistic Winds Derived 
Using "D" Value Approach with Observed 
Ballistic Winds (SIRS-B) 
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( 21 ) 



Ug = 152. 2 x 2 - 255. 8x 5 + 344* 7x 4 + 197.6x 6 



-695. 7x ? + 292 . 7x q + 14.74 



v B = -133.4y 2 - 121. Oy^ + 144. 4y 4 - 146. 0y 6 



+568. 8y ? - 129.6y 8 + 5595. 



( 22 ) 



As before, the x^ represent -j* and the y^ represent -j- . The 

results using these equations to calculate Ug and v.g are given in 
Table VIII, Again, the regression technique gave relatively poor 
results for the calculation of Ug in the independent sample. The 
high correlations for Ug in the dependent sample indicate that the 
low correlation for the independent sample may be due to data short- 
comings as previously indicated. Comparison with Table IV shows 
that the ±*e suits for v~ are better than for the January 3 IKS -A 
case. As was noted in the SIRS-A experiment, this method explained 
a somewhat lower percentage of the variance than did the "D" value 
approach. 



D. DIRECT APPROACH 

Direct correlations of the seven SIRS channel radiances with 
both Ug and v^ were calculated and the resulting regression 
equations were used as with the SIRS-A data. The results are 
tabulated in Table IX, Note that in this approach the sample sizes 
are much larger than in either of the other methods. 

This method also specifies Ug much less accurately than it does 
v B , The results for both components are not as good as those obtained 
with the SIRS-A data, (See Table V. ) 
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TABLE VIII 
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TABLE IX 



Comparison of Ballistic Winds Derived 
by Direct Approach with Rawinsonde-Derived 
Winds (SIRS-B) 



Dependent 

Sample 


n 


Corr 

Coeff 


Explained 

Variance 


Std 

Error 


Max 

Error 


u component 


130 


.288 


8.3% 


21kts 


-47kts 


v components 


130 


• 535 


28.6 


20 


-52 


Independent 

Sample 












u component 


64 


.302 


9.1 


17 


41 


v component' 


64 


.452 


20.4 


13 


30 



Because of the wide geographical distribution of the SIRS-B data 
the data sample is much more heterogeneous. The results of this 
test seem decisively insignificant. 

E. SUMMARY 

In neither the Ug nor Vg case were the results of the SIRS-B 
experiments as good as those obtained with the December SIRS-A data. 

For the v component the results were considered quite good considering 
the lack of window channel data and the wide spatial and temporal 
distribution of the data. Since all three approaches gave significantly 
poorer results for the u components, it is thought that a major portion 
of the difficulty must lie either in the nature of the u field or 
in the small size of the independent sample. 

When correlations between the rawinsonde-derived ballistic winds 
and winds geostrophically calculated from analyzed rawinsonde-derived 
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"D" values were computed, the results were mixed. For the entire 
test period (n = 45) "the correlation coefficient for u^ was .848 and 
for Vp was ,815» Making a similar test of correlation using only the 

J3 

data of Jan, 15-17 (n = 15) > however, the correlation coefficient for 
Ug was .812 and for Vg was .431* These anomalous results indicate 
the variability of the data in the period of the SIRS-B experiment 
and the need for more closely spaced SIRS data points in future 
experimentation. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Both the "D" value and gradient of radiance approaches were 
found to he capable of specifying the ballistic wind within the 
error limits set down by Finke [3] . Generally, the "D" value 
approach explained a higher percentage of the variance and gave 
lower standard errors. 

Contributing factors to the generally poor results for the u 
components upon independent data testing in the SIRS-B experiment 
may have been one or a combination of the following. The nature 
of the u-g field could be such that it is not well specified by a 
linear regression equation. It was noted that the regression 
equations usually failed to calculate any negative values of Ug. In 
only one of the two cases where negative values were calculated 
was the Ug obtained from rawinsondes also negative. Also, the 
previously exhibited poor temporal stability of the regression 
equations could be more pronounced for Ug in this data sample. 

The direct approach is not thought to be capable of claculating 
ballistic winds accurately enough to be useful. 



B. RECOMMENDATIONS 

It is recommended that experimentation be continued with a new 
SIRS-B data set having a more even data distribution around the 
northern hemisphere. Better data coverage would enable the use of 
shorter periods of time to derive the equations which would help 
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counteract the temporal stability problem* Additionally, rocketsonde 
data should be used in addition to rawinsonde data to more accurately 
define the motion and mass fields at upper levels. This, however, 
tends to limit data to the northern hemisphere. 

Experimentation with a long data record applied to homogeneous 
latitude should be conducted to determine the optimum time period 
for determining the regression equations. An optimum record length 
should be found which balances sample size and representativeness of 
the early portions of the record. 

In further SIRS-B experimentation, a more sophisticated method 
for calculation of gradients should be utilized. One method 
suggested by Kung [4] is to find the least squares approximation to 
the solution of a system of equations of the following form: 

= 4X; (25) 

Here Q, is any scalar variable, and Ax and Ay are the eastward 
and northward distances between the i XI pairing of the master data 
point and a slave. The unknown partial derivatives may be calculated 
if i 2 2 in each master-slave cluster. The greater the degree of 
overspecification, the better should be the least squares fit. 

Another method for obtaining gradients might be the description of 
the scalar surface as a polynomial by the of a regression equation 
in x and y. Derivatives could readily be calculated from the poly- 
nomial. This method would be time consuming since a separate 
regression equation would have to be generated for each set of points. 
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